Electronic structure of the strongly hybridized ferromagnet CeFe2 
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^ ^ . We report on results from high-energy spectroscopic measurements on CeFe2, a system of particular 

interest due to its anomalous ferromagnetism with an unusually low Curie temperature and small 
magnetization compared to the other rare earth-iron Laves phase compounds. Our experimental 
^ ■ results, obtained using core-level and valence-band photoemission, inverse photoemission and soft 

c 2 ^ ' x-ray absorption techniques, indicate very strong hybridization of the Ce 4/ states with the delocal- 

^3 , ized band states, mainly the Fe 3d states. In the interpretation and analysis of our measured spectra, 

C_H ■ we have made use of two different theoretical approaches: The first one is based on the Anderson 

impurity model, with surface contributions explicitly taken into account. The second method con- 
sists of band-structure calculations for bulk CeFe2 . The analysis based on the Anderson impurity 
model gives calculated spectra in good agreement with the whole range of measured spectra, and 
reveals that the Ce 4/ - Fe 3d hybridization is considerably reduced at the surface, resulting in even 
^ . stronger hybridization in the bulk than previously thought. The band-structure calculations are ab 

initio full-potential linear muffin-tin orbital calculations within the local-spin-density approximation 
of the density functional. The Ce 4/ electrons were treated as itinerant band electrons. Interest- 
i-pj ' ingly the Ce 4/ partial density of states obtained from the band-structure calculations also agree 

r^ ' well with the experimental spectra concerning both the 4/ peak position and the 4/ bandwidth, 

Q , if the surface effects are properly taken into account. In addition, results, notably the partial spin 

O ' magnetic moments, from the band-structure calculations are discussed in some detail and compared 

to experimental findings and earlier calculations. 
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PACS numbers; 71.27.-f-a, 75.75.Bb, 75.30.Mb, 79.60.Bm 

I. INTRODUCTION terns. In such systems, many unusual phenomena are 

typically observed, like for instance anomalously low sat- 

The 4/ states of rare-earth elements in solids usually uration magnetization and Curie temperature Tc (e.g., 

retain free-ionic properties with a well-defined integer CeFes), intermediate valence (e.g, SmS), heavy fermion 

occupation number. However, there are also rare-earth behavior (e.g., YbBiPt), or non-Fcrmi-liquid behavior 

compounds where the hybridization of the 4/ states with (e-g-, CeCug-^Au^). Even more surprisingly, simulta- 

extended band states is important, in which case they neous magnetic ordering and superconductivity has been 

may exhibit properties usually only found in actinide sys- observed (e.g., CeCuzSiz). The superconductivity is un- 



conventional, i.e., the order parameter suggests a d-wave 
superconducting state, as opposed to the conventional 
s-wave state. Finding a proper theoretical description 
of the 4/ states in these compounds remainS|-pne of the 
major problems in condensed matter physics.u 

As for the photocmission spectroscopy (PES) studies 
of Ce compounds, it is widely believed that the spectra 
are well described by the single-impurity Anderson model 
(SIAM).H Recently, however, it has been argued that sys- 
tems, in which the Ce 4/ states hybridize strongly with 
the other valence electrons, calculations based on den- 
sity functional theory (DFT) may give an equally good, 
or even better description of the photo-emission spec- 
tra than the SIAM analysis, provided that surface^ffects 
are properly taken into account in the analysis .a How- 
ever, one should bear in mind that calculations based on 
DFT are not strictly applicable for excited state prop- 
erties, instead the ground state properties, such as the 
magnetic moments of the ground state, which are typ- 
ically the focus of these calculations. Nevertheless, the 
electronic structure given from such calculations are of- 
ten compared with photocmission data and good agree- 
ment between experiment and calculations is frequently 
observed. In the limit of complete screening of the excited 
state, one would expect ground-state density functional 
calculations to be able to describe the spectra well. 

CeFe2 is thought to belong to a class of strongly hy- 
bridized systems. This compound shows ferromagnetism 
below Tc = 230 K with a saturation magnetization of 
2.30 fiB/i-'^- Above Tc, the magnetic susceptibihty fol- 
lows the. Curie- Weiss law with an effective moment of 7.4 
/is/f.u.Q If one compares CeFe2 with the other i?Fe2 com- 
pounds (R: rare earth elements) , a number of anomalies 
in its physical properties can be observed. The lattice 
constant is much smaller than an interpolation using the 
lattice constants of the other it!Fe2 systems would sug- 
gest. Its Curie temperature is anomalously low: the other 
_RFe2 compounds have Curie temperatures ranging from 
596 K to 796 K.El The saturation magnetization is unusu- 
ally low compared to the other i?Fe2 compounds (2^3 
and 2.90/is/f.u. for LuFe2 and YFe2, respectivelyD). 
Moreover, even if only a small fraction of the Fe atoms 
are substituted for Al, the ferromagnetic ordering is di 
stroyed, and the system becomes anti-ferromagnetic.l 
In fact, even in pure CeFe2, recent neutron scattering 
experiments have revealed strong competition between 
the ferromaguetic ground state and an antiferromagnetic 
ground statcQ Together, these facts suggest that the Ce 
4/ states in CeFe2 hybridize strongly with the other va- 
lence electrons, notably the Fe 3d valence states. This 
hypothesis is further supported by_the X-ray absorption 
(XAS) experiments by Croft et al.S 

In this paper, we present high-energy spectroscopic 
results on CeFe2 including core-level x-ray photocmis- 
sion (XPS), XAS, Ce 3<i-4/ and 4d-4/ resonant PES, 
bremsstrahlung isochromat (BIS) and high-resolution ul- 
traviolet photocmission spectroscopy (UPS) in order to 
elucidate the electronic structure of this system. In the 



case of Ce compounds with strongly hybridized 4/ states, 
it has been pointed out that surface effects are. extremely 
important in the interpretation of the spectra. Bll^l There- 
fore, we have attempted to differentiate the electronic 
structure of bulk and that of surface for CeFe2 in the 
analysis o f the s pectra. As will be further elaborated on 
in section [II C of this paper, electronic-structure calcu- 



lations with the Ce 4/ states treated as valence s 
give a good description of the magnetism in CeFe2]i3'ta 
It is of course highly interesting to assess the applica- 
bility of the same theory in describing also the photoc- 
mission spectra of CeFe2, even though as noted these 
calculations are not strictly applicable for|-epccited state 
properties. Very recently, Sekiyama et a^Jlil reported a 
high-resolution 3rf-4/ resonant photocmission study of 
the strongly hybridized system CeRu2 and found that the 
Ce 4/ spectra can be explained by band theory. In the 
following, we first attempt to describe the spectra in the 
framework of the SIAM and obtain a set of SIAM param- 
eters. In the analysis, surface effects on the spectra are 
explicitly taken into account. Next, the bulk component 
of the valence-band spectra is compared with the density 
of states calculated using density functional theory. All 
DFT results presented here have been calculated within 
the local spin-density approximation (LSDA) . The use of 
more recently developed generalized gradient functionals 
would not, however, alter any of our conclusions. 



II. METHODS 

A. Experiment 

Polycrystalline samples of CeFe2 were prepared by arc- 
melting the pure constituent materials. Subsequently, 
the samples were annealed at 750 °C for a week to ob- 
tain single phase samples. Magnetization measurements 
yielded the same Tc as in the literature. The XPS spec- 
tra were taken with Mg Ka radiation {hi/ — 1253.6 eV) 
using a double-pass cylindrical-mirror analyzer, and the 
BIS spectra were obtained using a Pierce-type electron- 
gun and a quartz crystal monochromator which was set 
at hv = 1486.6 eV. The Ce 4d-4/ resonant PES mea- 
surements were done at beam-line BL-2 of SOR-RING, 
Institute for Solid State Physics, University of Tokyo. 
The Ce 3d-4/ resonant PES and Ce 3d XAS data were 
taken at beam-line BL-2B of Photon Factory, High En- 
ergy Accelerator Research Organization. Photoelectrons 
were collected using a double-pass cylindrical-mirror ana- 
lyzer in the resonant PES measurements. The XAS spec- 
tra were obtained by measuring the total electron yield 
using an electron multiplier placed near the sample. All 
measurements were done in the range 50-80 K, i.e., be- 
low the Curie temperature. In the case of the Ce Ad-Af 
resonant PES, additional measurements at room temper- 
ature, i.e, above the Curie temperature, were performed. 
The total energy resolution was '^1.0 eV for XPS and 



BIS, ~0.5 eV for Ce 4d-4/ resonant PES, ~0.5 eV for 
XAS and -1.0 eV for Ce 3d-4/ resonant PES. The high- 
resolution UPS measurements were done around 17 K 
using a hemispherical analyzer and the He I (hv — 21.2 
eV) and He II {hi' = 40.8 eV) resonance lines. The energy 
resolution was —25 meV for both photon energies. The 
binding energies were calibrated using Au evaporated on 
the samples. For XAS and Ce 3(i-4/ resonant PES, the 
photon energies were calibrated using the Cu 2p edge of 
Cu metal and the Co 2p peak of LaCoOs. Clean surfaces 
were obtained by scraping the sample repeatedly, while 
maintaining the sample under ultra-high vacuum, with a 
diamond file prior to each measurement. Cleanliness of 
the surfaces was checked by the absence of O Is and C Is 
XPS signals from contaminants in the case of the XPS, 
XAS, BIS and Ce 3d-4/ resonant PES measurements. In 
the case of the high-resolution UPS and Ce Ad-Af reso- 
nant PES measurements, cleanliness was checked by the 
absence of a O 2p feature which appears around 6 eV 
below the Fermi level [Ep)- 



B. Single-impurity Anderson Model 

The SIAM calculations were made based on the vari- 
ational l/A^/-expansionr4jjiethod developed by Gunnars- 
son and SchonhammerJiJ Here, we performed the cal- 
culations to the lowest order in l/Nf, where Nf is the 
degeneracy of the Ce 4/ level and was taken to be 14. 
The /^ configuration was also included in the calcula- 
tion. The energy dependence of the hybridization ma- 
trix elements was taken from the off-resonant spectra, 
which approximately represent the Fe 3d partial density 
of states. The configuration dependence of the hybridiza- 
tion strength was also taken into account, and was chosen 
to be the sanie as that obtained for a-Ce by Gunnarsson 
and Jepsen.Ej In the calculations, we divided the ba»d 
continuum into discrete levels following Kotani et alxl 
We further assumed that each spectrum was a superposi- 
tion of two components which represent bulk and surface 
spectra. The weight of each component was treated as 
fitting parameters within a range consistent with the uni- 
versal curve for the mean free path of photoelectrons.Efl 



C. Band-Structure Calculation 

In the band-structure calculations presented here, we 
have used the fuU-natential linear muffin-tin orbital 
method (EP-LMTO).EI In this method, the Kohn-Sham 
equationaH3 are solved for a general potential without any 
shape approximation. The local (spin) density approx- 
imation (LSDA) in the Hedin-Lundqvist parameteriza- 
tiorO was used for the density functional. 

In the FP-LMTO method, space is divided into non^ 
overlapping spheres, so-called called muffin-tin spheres,LJ 
surrounding each atomic site, and an interstitial region. 



The basis functions used are energy independent Bloch 
functions, whose construction is different in the spheres 
and in the interstitial. 

A basis function in the interstitial is defined by the 
Bloch function of solutions to the spherical Helmholtz 
equation with nonzero kinetic energy k^, or a linear com- 
bination of such solutions for different kinetic energies. 
The Fourier representation of this basis function is taken 
from the Fourier series of a function matching the ba- 
sis in the interstitial region but not inside the spheres, a 
so-called pseudo-wave function, whose exact shape inside 
the muffin-tin sphere is of no importance for the final so- 
lution as long as it is continuous and differentiable at the 
sphere boundary and matches the true basis function in 
the interstitial. 

Inside the spheres, where the charge density varies 
rapidly, the basis functions are Bloch functions of nu- 
merical radial functions times spherical harmonics. The 
radial part of a basis function is constructed from the 
numerical solutions (pL {Eu , r) of the radial Schrodinger 
equation in a spherical potential at the fixed energy E^, 
and their energy derivatives 4'L{Ei,,r). Here, the index 
L stands for a collection of quantum numbers: the prin- 
cipal quantum number n, the orbital quantum number /, 
the magnetic quantum number m, and the kinetic energy 

The treatment of the entire basis set within one sin- 
gle energy panel allows all states, including the semi-core 
states, to hybridize fully with each other. Our method 
is linear, i.e., the basis functions are constructed by ex- 
panding around fixed energies E,j. The expressions for 
the crystal wave functions in the muffin-tin spheres are 
matched to the interstitial crystal wave function at the 
sphere boundaries so that the total crystal wave function 
becomes continuous and differentiable in all space. In the 
present calculation, the expansion in spherical harmonics 
was taken up to Z = 6. For Ce, the 6s, 5p, 6p, 5(i, and 4/ 
orbitals were included in the basis set, with 5p as semi 
core. For Fe, we included the 4s, Ap and 3d orbitals, i.e., 
no semi-core state was used for Fe. Four K^-values were 
used in the calculation: -0.6 Ry and -0.1 Ry for the va- 
lence states, and -1.5 Ry and -1.0 Ry for the semi-core 
Ce 5p states, all with respect to the muffin-tin zero. 

Reciprocal space was sampled with what would corre- 
spond to 1331 k-points in the full Brillouin zone (BZ) 
using special k-point sampling methodsHj The non- 
overlapping muffin-tin spheres were chosen as 21% and 
17% of the unit-cell lattice constant for Ce and Fe, respec- 
tively. With this choice, 36% or the unit cell volume is in 
the interstitial region and the closest muffin-tin spheres 
are 3% from touching. 

The experimental lattice constant was used in the 
calculations. Furthermore, the calculations were spin- 
polarized but the spin-orbit interaction was not included. 
This latter approximation will be commented on further 
in conjunction with presenting and discussing the results 
from the band-structure calculation. 
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FIG. 1. Core-level spectra of CeFe2. (a) Ce 3d core-level 
XAS spectra, (b) Ce 3d XPS spectra taken at hi/ = 1253.6 
eV. 



III. RESULTS 

A. Experiment 

Fig. ^ shows the Ce 3d core-level XPS and XAS spec- 
tra. The XPS line-shape is a typical one for a strongly 
hybridized Ce compound, consisting of three peaks which 
correspond to the StfAf, 3^4^ and 3^^4/2 final 6dat£s 
in each of the j—3/2 and 5/2 spin-orbit componcntsLllt3. 
In the XAS spectrum, the main peaks are due to the 
3d^4/^ final-state multiplet and the satellite structures 
^--^5 eV above the main peaks are due to the 3d^4/^ final 
states. The rather distinct 4/'^ peaks in the XPS spec- 
trum and the 4/^ structures in the XAS spectra, together 
with the obscured 3(PAf^ final-state multiplet structures 



of the main XAS peaks, indicate strong hybridization of 
the 4/ states with the valence band in this system. The 
XPS spectrum reflects the surface electronic structure be- 
cause of the rather low kinetic energies of photoelectrons 
from the Ce 3d core level. A detailed analysis of this is 
given below. 
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FIG. 2. Valence-band PES and BIS spectra of CeFe2. 

hv = 121 and 114 eV (881 and 875 eV) correspond to Ce 4d-4/ 

(3d-4/) on- and off-resonance, respectively. Solid curves show 

the difference spectra, which represent the Ce 4/ component. 

The results of valence-band PES and BIS are shown in 
Fig. ^. The on- and off-resonance occurs, respectively, at 
hv ^ 121 and 114 eV in the Ce 4d-4/ resonant PES and 
at hv = 881.1 and 875.4 eV in the Ce 3d-4/ resonant PES. 
Identical spectra in the present resolution have been ob- 
tained for Ce 4d-4/ resonant PES at room temperature, 
which is above Tc (not shown). We have obtained the 
Ce 4/ spectra by subtracting the off-resonance spectra 
from the on-resonance spectra as shown by solid curves 
in Fig. g. As seen from this figure, there is a large dif- 
ference between the Ce 4/ spectra obtained from the Ce 
4d-4/ and Ce 3d-4/ resonant PES. While the former has 
a double-peak structure in the vicinity of -2 eV and near 
Ep, the latter is dominated by a single peak near Ep, 
implying stronger hybridization in the latter. This can 
be attributed to the difference in the surface sensitivity 
of the two spectra due to the different kinetic energies of 
photoelectrons. This also indicates that the Ce valency 
at the surface is closer to trivalent than it is in the bulk. 

In Fig. H we also show the XPS spectrum of the va- 
lence band taken with Mg Ka radiation. Owing to 
the higher kinetic energies of photoelectrons, this spec- 
trum is considered to be more bulk sensitive than the 
above PES spectra. Considering the photoionization 
cross-sectionsjHj the valence-band XPS spectrum should 
mainly reflect the Fe 3d partial density of states (DOS) 
with significant contributions from Ce 4/ and Ce hd. As 
seen in the figure, the XPS spectrum shows a line-shape 
similar to the off-resonance spectra of Ce Ad-Af and Ce 



3d-4/ resonant PES. 

It is expected, that the BIS spectrum should also re- 
flect the bulk electronic structure rather well. There is 
a peak near Ep and a broader feature at ~6 eV. They 
originate mainly from the Ce 4/ states, although there 
are contributions from the Ce Srf and Fe 3d states, too. 
The peak near Ep and the structure around '^6 eV coiu 
respond to the 4/^ and if^ final states, respectively.El 
The broad line-shape of the gtructure at '-^6 eV is due to 
the 4/^ final-state multiplet.D The strong intensity of the 
peak near Ep again indicates strong hybridization of the 
Ce 4/ states with the valence states. 



of the mean free path of photoelectrons.EJ In the near 
Ep region, structures just below Ep and at ~ -0.3 eV 
are somewhat enhanced in the 40.8 eV spectrum. These 
structures originate from the Ce 4/ states and correspond 
to the tail of the Kondo resonance (possibly with unre- 
solved fine structures due to crystal-field splitting) and 
the spin-orbit side band, respectively. These structures 
are also expected to be dominated by surface contribu- 
tions. 



B. Single- Impurity Anderson Model 
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FIG. 3. High-resolution UPS spectra of CeFe2. Inset 
shows an enlarged view near Ep- 

Fig. H shows high-resolution UPS data. In this photon 
energy range, the cross sections of the Fe 3d, Ce AI and 
Ce 5d states varies rapidly with photon energy.Ej The 
relative cross-section of Ce 4/ to the other orbitals in- 
creases when going from hiy =21.2 to 40.8 eV while that 
of Ce 5d rapidly decreases. Therefore the structure at -(2 
~ 3) eV which appears only in the hv = 40.8 eV spec- 
trum, originates from the Ce 4/ states, and corresponds 
to one of the double peaks in the Ce 4/ spectrum ob- 
tained by the Ce 4d-4/ resonant PES. This observation 
is also consistent with the fact that the 40.8 eV spectrum 
is surface sensitive, according to the "universal curve" 



The SIAM parameters obtained in our calculation are 
listed in Table |. Here, e/ is the position of the bare 
4/ level (4/^ -^ 4/° ionization level) relative to Ep, Ug 
is the Af-Af on-site Coulomb energy, C//c is the 4/-core- 
level Coulomb energy and V is the Ce 4/-valence-band 
hybridization strength, in accordance with the definitions 
in Ref. O. Using those parameters, the 4/ occupation 
number Uf has been calculated and listed in the last col- 
umn of Table H. Since we have fitted many different types 
of spectra using the SIAM shown in Fig. ^, many con- 
straints have lead to a rather unique set of SIAM param- 
eters. 

In Fig. ^, comparison is made between the SIAM cal- 
culations and the experimental spectra. As seen, we ob- 
tain good overall agreement with all experimental spec- 
tra. The main discrepancy between the SIAM results 
and experimental spectra is found in the BIS spectrum 
(Fig. He), where the position of the calculated /^ peak is 
about 0.5 eV lower than in the experimental spectrum. 
Noticeable from this figure is also the large difference be- 
tween the bulk and surface spectra obtained through the 
SIAM analysis. For instance, in the Ce 4/ spectrum ob- 
tained from Ce Ad-Af resonant PES (Fig. |b), which is 
basically a double-peak structure, the relative strength of 
the two peaks is very different between the bulk spectrum 
and the surface spectrum. Similar differences in directly 
measured bulk- and sucface spectra have previously been 
reported for Ce-metal,E2l and caH-ie explained as due to 
larger hybridization in the bulk.Ej 

Tablel shows that, apparently, the Ce atoms belonging 
to the surface are as good as completely trivalent, with 
Uf ~ 1.0, whereas in the bulk, the 4/ states are strongly 
hybridized, having the significantly lower occupation of 
0.78. 



TABLE I. 

Ufc are given 


SIAM parameters 
in units of eV. 


for CeFe2. 


£/, Uff, V, and 




^/ Uff 


V 


Ufc Uf 


surface 
bulk 


-1.8 6.4 
-0.8 6.4 


0.23 
0.41 


9.7 1.0 
9.7 0.78 



However, there are noticeable amplitudes of the /" 
and /^ configurations also at the surface, indicating that 
also here, some hybridization between the 4/ and valence 
states is taking place. 
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FIG. 4. Comparison of the single-impurity Anderson 

model calculation with the experimental spectra of CeFe2. 

(a) Ce 4/ spectrum obtained from Ce 4d-4/ resonant PES. 

(b) Ce 4/ spectrum obtained from Ce 3d-4/ resonant PES. 

(c) BIS spectrum, (d) Ce 3d core-level XAS spectrum, (e) Ce 
3d core-level XPS spectrum. In each panel, dots show experi- 
mental spectrum, solid curve shows the calculated spectrum, 
and dotted and dash-dotted curves show the calculated sur- 
face and bulk components, respectively. 



C. Band-Structure Calculation 



Experimentally, the partial moments in CeFe2 
have been studied using several different experimental 
methods: polarized neutronsJ£3 Compton scatteringEEl 
and, verjj— recently, x-ray magnetic circular dichroism 
(XMCD).E3 In all experiments, an anti-parallel coupling 
of the Ce and Fe moments is found. This coupling is 
also reproduced .in our calculation, jEs well as in earlier 
calculations .E3E3 As is well known Jl3 this anti-parallel 
coupling of the moments is a strong indication that the 
Ce 4/ states in CeFe2 are delocalized. This can easily 



be understood from the following argumentation. If the 
4/ electrons are localized, the 4/ spin moment would be 
dictated by the polarization of the spd-electrons of the 
Ce atom, that via hybridization effects are known to be 
anti-parallel to the 3d moment of the Fe atom. Hence 
the syin moments of the Ce atom and the Fe atom are 
always anti-parallel, both in the localized and delocalized 
case. For localized 4/ electrons, the 4/ spin moment is 
accompanied by an orbital moment (larger than the spin 
moment) that (via Hund's third rule) is anti-parallel to 
the Ce spin moment. Hence, for localized 4/ electrons 
the total (spin-f orbital) Ce-Fe coupling is ferromagnetic, 
whereas if the Ce 4/ orbital moment is quenched, due to 
band formation, the coupling is anti-parallel. 

In Fig. 1^, the spin-resolved partial DOS for the Ce 4/, 
Ce hd and Fe Sd states are shown. The first and third 
panels show the majority spin channel for Ce and Fe, 
respectively, and the second and fourth panels show the 
minority spin channel. Comparing the DOS for Ce and 
Fe, we see that the Ce 4/ and Fe 3d states have opposite 
spin polarization. Furthermore, the Ce 5d band width is 
seen to be much larger than that of the Ce 4/ and Fe 3d 
states, with the magnitude of the Ce 5d DOS roughly an 
order of magnitude smaller than that of the Ce 4/ and 
Fe 3d states. 

Fig. y shows the band-structure of spin-polarized 
CeFe2 along high-symmetry directions in the Brillouin 
zone. The flat bands clustered just above the Fermi level 
are predominantly of 4/-character. In the region from 
^2 eV to ~I0 eV, and around -5 eV with respect to the 
Fermi level, the spin splitting of the bands is clearly vis- 
ible. 

Our calculated total spin magnetic moment amounts to 
2.48 /iB per formula unit (experimental saturation mag- 
netization: 2.30 /iB, as stated earlier in this paper), with 
the main contributions being the Ce 4/ moment -0.54 
/iB, the Ce 5d moment -0.23 /iB, and the Fe 3d moment 
1.75 /iB- The partial occupation numbers summed over 
spin, within the muffin-tin spheres, are 1.07 for Ce 4/, 
1.28 for Ce 5d, and 6.18 for Fe 3d. Note that the partial 
spin magnetic moments are calculated using the partial 
occupation numbers inside the muffin-tin spheres, which 
is somewhat arbitrary. The total spin moment is on the 
contrary, of course, well defined. An obvious point, which 
seems to have been overlooked so far, is that not only in 
band-structure calculations, but also experimentally, the 
division of space between individual atomic species in a 
compound is in fact not unique nor even well defined. It 
is reasonable to believe that different experimental proce- 
dures differ in their "volume of sensitivity" around each 
atom, and thus effectively correspond to different ways 
of dividing up the total space in the compound between 
the atoms. This could be one reason why different exper- 
imental techniques find quite different values for the par- 
tial magnetic moments, and also why, in order to find the 
total Ce moment from experimental results, assumptions 
have to be made regarding the ratio of the number o£,5d 
and 4/ electrons contributing to the magnetization.EZi 
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FIG. 5. Spin-resolved partial DOS for the Ce 4/, Ce 5d 
and Fe 3d states. The Fermi energy is taken as the energy 
zero. The first and third panels show the majority spin chan- 
nel, and the second and forth panels, the minority spin chan- 
nel. In order to enhance visibility, the magnitude of the Ce 
5d DOS (dashed line) has been multiplied by a factor of 10. 
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FIG. 6. Energy bands for CeFe2 along high-symmetry 
directions. The Fermi energy is taken as the energy zero. 

An analysis of experimental data along this direction 
of thought might help resolve jControversies regarding the 
electronic structure of CeFe2 .EJ 

A calculation of the moments-iucluding spin-orbit cou- 
pling and orbital polarization,^ using the FP-LMT-Q 
method, has been performed earlier by Trygg et aZ..t3 
The difference between the presently reported fftin mo- 
ments and the ones reported by Trygg et aL,t3 which 
include spin-orbit coupling, is very small, around 1%. 
Thus, the effect of including spin-orbit coupling is shown 
to have only a very minor effect on the magnitude of 
the spin moment. Calculations by Eriksson et alJt3 us- 
ing the atomic sphere approximation (ASA) give some- 
what different values for the magnetic moments than the 
present method, in which no such geometrical approxima- 
tion regarding the form of the potential, wave functions 
or charge density is made. As demonstrated in Ref. nq, 
the spin density in CeFe2 is highly non-spherical, which 
may well be the reason for the differences in results from 
full-potential and ASA calculations. To summarize, the 
arguments presented above justify our present calcula- 
tional approach, i.e, using a full-potential method, but 
neglecting spin-orbit coupling. 

Concerning the absolute magnitudes of the individual 
moments, the discrepancies between different experimen- 
tal approaches can be quite large, for instance, the Ce 4/ 
spin moment is measured to be -0.37 fis with XMCD, 
whereas polarized neutrons find the corresponding mo- 
ment to be only about a fourth as large: -0.10 ^b- As 
already touched upon above, one reason for these discrep- 
ancies between different experimental techniques may 
well be that they differ in the way the space in the com- 
pound is effectively divided up between the atoms. With 
this in mind, the magnitudes of our calculated moments 
must be said to be in satisfactory agreement with exper- 



imental findings, although the overall trend appears to 
be that the calculations overestimate the moment mag- 
nitudes. 
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FIG. 7. Comparison of the DFT DOS with the experi- 
mental spectra, dot-dashed curves show orbital components. 
In the PES part, the Ce 4/ spectrum obtained by Ce 3d-4/ 
resonant PES is shown. 

In Fig. M, we compare the valence-band Ce 4/ PES and 
BIS spectra with the DFT DOS. As mentioned earlier, 
the Ce 3d-4/ resonant PES spectrum and the BIS spec- 
trum are rather bulk sensitive, and thus it is relevant 
to compare these spectra with bulk DFT calculations. 
In the PES part of the spectra, comparison is made be- 
tween the experimentally obtained Ce 4/ spectra and the 
Ce 4/-projected DFT DOS. As for the BIS part, the Ce 
4/, 5(i and Fe 3d partial DOS have been added takiug 
account of the atomic photoionization cross-sections .E3 
Agreement between experiment and theory is satisfac- 
tory almost to the same extent as in the SIAM calcu- 
lation. In the DFT DOS, the structure around 6 eV in 
the BIS spectrum is of course not reproduced, since this 
structure corresponds to the 4/^ final state, and thus is 
a purely excited-state property of the system. Further- 
more, in the DFT calculation, the energy of the near Ep 
peak in the Ce 4/ PES spectrum and also in the BIS spec- 
trum is slightly higher than in the experimental spectra. 
Also, the intensity on the higher binding energy side of 
the PES spectrum is underestimated in the calculation. 
However, one should note that the relative intensities de- 
pend also on the transition matrix elements, which are 
not included in the DOS curves. 

Fig. 11 shows the valence-band XPS, Ce Sd-Af and Ad- 
4/ off-resonance and UPS (He II) spectra. The spec- 
tral weight comes primarily from the Fe 3d states, and 
thus we compare these spectra with the Fe 3d partial 
DFT DOS. Although there are differences in the surface 
sensitivity and in the contributions from other orbitals, 
all the experimental spectra have similar band widths 
and line shapes. In comparison with the DFT DOS, al- 
though overall features are well reproduced, the experi- 
mental spectra have larger spectral weight near Ep than 
the calculated DOS. 
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IV. DISCUSSION AND CONCLUSIONS 

Generally, for any material, the surface electronic 
structure can differ substantially from the bulk one. For 
valence fluctuating systems like CeFe2, it is very likely 
that the electronic structure of the Ce atoms close to the 
surface is not the same as that of the bulk Ce atoms. 
Therefore, in order to study the bulk electronic structure 
of valence fluctuating Ce compounds by means of high- 
energy spectroscopic methods, it is essential to take into 
account the effects of the sample surface when interpret- 
ing the spectra. In the present work, by assuming that 
the spectra are superpositions of the surface and bulk 
components, we have shown that all measured spectra of 



CeFe2 are fairly well reproduced by the SIAM calcula- 
tions using the same set of parameters. At the surface, 
Ce is found to be nearly trivalent. The bulk set of pa- 
rameters places CeFe2 in the strongly intermediate- valent 
regime, giving a 4/-occupation number Uf as small as 
0.78 (with considerable amplitude of the p configura- 
tion both in the surface and bulk). This means that the 
Ce 4/ states are strongly hybridized with the Fe 3d states 
in the bulk and that the states around Ep have a large 
amount of / character. Also, the differences between 
our SIAM-derived bulk- and surface spectra for CeFe2 
are similar to the differences between directly measured 
bulk- and surface spectra of Ce-metal, a difference which 
can be explained as due to larger 4/ hybridization in the 
bulk than at the surface. 

Apart from the SIAM analysis, a number of features 
of the measured spectra force us to draw the same con- 
clusion regarding the nature of the 4/ states in the bulk 
and at the surface, notably the large difference between 
the Ce 4/ spectra obtained from the Ce 3d — 4/ and Ce 
4d — Af resonant PES, and the strong intensity of the 
peak near Ep in the BIS spectrum. 

As a result of the strong hybridization of the Ce 4/ 
states in the bulk, the "/°" final state feature in the Ce 
4/ spectrum deduced from the Ce 3d-Af resonant PES 
has a very weak intensity in contrast to the Ce Ad-Af 
resonant PES spectrum, where the /''-final-state feature 
leads to the well-known double-peak structure. In such a 
case, i.e, where there is strong screening of excitations, it 
is expected that the bulk 4/ spectrum can be interpreted 
in terms of a one-electron picture, and thus the DFT DOS 
should compare well with the 4/ spectrum. This is also 
seen to be the case. Comparison of the 4/ spectra with 
the DOS (Fig. I) shows that the Ce 4/ partial DFT DOS 
describes the valence-band Ce 4/ spectra well, except for 
the /^ structure in the BIS spectrum of course, since 
the /^-peak is due to incomplete screening. This fact 
poses the question of how the "/^" final state (which 
is commonly referred to as the "Kondo peak") of the 
SIAM picture and "the 4/ band" in the band picture are 
related to each other, since according to DFT, this peak 
is a one-electron feature, and in the SIAM, this is due to 
a many-body effect. 

We also draw the conclusion that due to the strong 
hybridization in the bulk, the spin-orbit side band seen 
in the high-resolution UPS spectra (Fig. ^) must have 
surface origin since it is known that the spectral weight 
of the spin-orbit side bpid is strongly reduced when the 
hybridization is strong.113 

We now turn to a more detailed comparison of how 
well the DFT calculations and the SIAM analysis perform 
for the different spectra. Regarding the position of the 
near Ep peak in the BIS spectrum, the DFT calculation 
predicts a higher energy than experiment while the SIAM 
calculation predicts a lower energy than experiment. The 
intensity on the higher binding energy side of the near 
Ep peak in the PES spectrum is underestimated in the 
DFT calculation (although strictly speaking, intensities 



cannot be expected to be reproduced with a DOS, since 
the transition matrix elements are neglected) while it is 
overestimated in the SIAM calculation. The position of 
the near Ep peak in the PES spectra is calculated to be 
too close to Ep both in the SIAM and DFT calculations. 
As for the Fe 3d component, the DOS calculated using 
DFT does not give a good account of the intensity in 
the experimental spectra near Ep. This may be due to 
that the 4/ transition matrix elements are large close 
to the Fermi level compared to the 3d transition matrix 
elements. 

All in all, the above discussion amounts to that the Ce 
4/ electrons in the bulk hybridize strongly with the Fe 
3d electrons. This conclusion agrees perfectly with the 
experimentally observed anti-parallel coupling of the Ce 
and Fe moments, which is also reproduced in the DFT 
calculation. In the SIAM, the 4/ electron is assumed 
to be localized, which indirectly implies that a parallel 
coupling of the Ce and Fc moments is expected. 

Finally, we wish to mention some sources of error in the 
present work. Our measurements were done on scraped 
surfaces, which might make the surface rough, and there- 
fore ill-defined. Furthermore, the precise values of photo- 
electron mean free paths are difficult to estimate, which 
naturally also has the effect of making the border between 
"bulk" and "surface" somewhat ill-defined. Our SIAM 
is not the most elaborate one, for instance we assume 
a degeneracy of 14 of the 4/ level, thereby neglecting 
spin-orbit coupling and anisotropic hybridization effects, 
which leads to a crystal-field splitting. Furthermore, as in 
all DFT calculations, the functional used treats electron 
correlation only to a limited extent, i.e., it is not meaning- 
ful to expect perfect agreement between the DFT results 
and experiment. 
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